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Abstract 

A ground state energy variational calculation of anyon gas with Hamiltonian included the in- 
teraction of spins of particles with anyon vector potential induced, i.e. statistical, magnetic field 
exhibits exact cancelation of terms connected with fractional statistics. This leads to bosoniza- 
tion of anyons due to coupling of their spins with statistical magnetic field. We presume that at 
the dense gas fluctuations of effective spins destroy the coupling and bosons become anyons. At 
the assumption that pseudogap (PG) boundary is temperature independent and when anyons are 
fermions we use this model to interpret experimental phase diagrams of Tallon and Loram hole and 
electron doped High- Tc superconductors below PG energy Eg and find the qualitative and quanti- 
tative agreement. We do the hypothesis that phase transition (PT) of bosons into Bose-Einstein 
condensate is not of second order, but of first order, close to second one, PG regime is meta stable 
phase of bosons, and Eg = is the critical point of this PT. Bosons undergo PT into fermions 
on PG boundary. Described in the literature non-Fermi quasi-particles might be related to bosons 
with effective spins. 
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I. INTRODUCTION 



In 2D systems the concept of anyons provides us with unique opportunity to introduce 
antisymmetric property of fermion wave function into Hamiltonian. It is assumed that 
particles are to be spinless. The background of physics of anyons, the Aharonov-Bohm 
effect, itself is interesting phenomenon. Related experimental investigation |l| has been 
recently performed. Anyons enable to explicitly investigate relationship of spins of particles 
and statistics. From standard courses of relativistic quantum theory (see, for example, 
Ref. 2|]) it is well known that particles with integer number oi h/2 spins should obey a 
Fermi statistics. It seems, this connection between spins and statistics is strong. In this 
work, however, we show that the gas of 2D anyons or fermions with spin h/2 can undergo 
unexpected changing of statistics and this occurs due to Zeeman interaction of spins with 
statistical magnetic field 3j produced by vector potential of anyons. The calculation of 
ground state energy expectation value will be carried out in the framework of variational 
approach with cut-off parameter regularization P, 0] , which we developed recently. 



II. THEORY AND MATHEMATICAL TREATMENT 



The Hamiltonian is 



H 



1 ^ 

— Y 

2M ^ 



k=l 



1 ^ 

4e 



k=l 



N 



(1) 



for gas of anyons with mass M and charge e, confined by 2D parabolic well, interacting 
through Coulomb repulsion potential, in the presence of uniform positive background 6| 
V{fj?). Here and pk represent the position and momentum operators of the A;th anyon in 
two space dimensions. 



N 



(2) 



is the anyon gauge vector potential [71], rfcj = r/c — r^-, and is the unit vector normal to the 
2D plane. The factor v determines the fractional statistics of anyon: z/ = (bosons) and 
V = \ (fermions). 
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In the variational scheme jj] we minimize the expression 

^ ^ / ^um^AR) dR 

J ^*j,{R)^t{R) dR ' 

Here R = {fi....fN} is the configuration space of the anyons. When energies are expressed 
in units of hLjQ = h?/ [ML"^) and lengths in units of L the normalized trial wave function in 
the bosonic representation of anyons reads 

N 



k=l ^ ^ 

Here a is variational parameter. The harmonic potential regularization with tending 
number of particles to infinity yields the ground state energy of infinite Coulomb anyon 
gas. 

Now we introduce in the Hamiltonian the term 

N 

M 



k=l 



with statistical magnetic field 



which can be derived if calculates 6^ = V x y4,^(rfc) by using Eq. ([2]). The sign in Eq. ([5]) is 
taken for electrons with charge e = — |e|. For holes, with charge e = |e|, we need to change 
a sign for z/ in Eqs. and then Eq. as also the expectation value for energy, Eq. 
([9]), (see below), will retain the sign. 

For Sz = h/2 and if we take into account that length unit is L, then (5*^^^(r) should be 
replaced by 5^'^\r)/L'^, hence, 

k=l k,j{k^j) 

The calculation of expectation value Eq. ([3]) when Hamiltonian is Eq. ([7]) and wave 
function ^>t{R) is Eq. (jlj) gives (we omit a factor h^/{ML^)) 

-Tiv I ^t{R) S^^\fk - rj) ^t{R) dR 
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The total expectation value for energy, Eq. ([3]), including all terms of Eq. (jT]), is 

E=^ + ^ + NM a^'^ ^ 1 , (9) 

where the term for M. is responsible for Coulomb interaction (see Refs. 0, S]). 

We did in a cut-off parameter regularization and found M = 1 + i/(A^ — 1). This expres- 
sion for A/" successfully describes energy of confined anyons with and without Coulomb 
interaction as well as one of infinite anyon gas [5] with Coulomb interaction. Substituting 
A/" into energy E, Eq. ([9]), we see the exact cancelation of terms with u factors. This result 
for energy can be obtained if we put z/ = 0, i.e., for case of bosons. As the energy of bosons 
is lower than one for anyons and fermions, there appears a coupling of spin with statistical 
magnetic field for every particle or bosonization of 2D anyons and fermions. 

One can assume the fluctuations of spins coupled to magnetic field. Therefore, bosons 
with effective spins might look like as Fermi particles. However, fermions with different spins 
are independent . Thus, the spins of bosons interact with each other and do not interact 
with spins of another fermions if they exist in the system. We introduce a some correlation 
length, inside of which spins of bosons interact with each other. For temperature T = 
we denote it C,o. The increase of fluctuations destroys the coupling, and bosons become the 
anyons or fermions. This occurs when the gain in the energy due to fluctuations of spins of 
bosons is equal to energy difference between the anyon (or Fermi) and Bose ground states. 

The interaction of spins of bosons we bring in the form 

N 

e-^°/«° . (10) 

fe=i 

Here it was introduced a factor e~^°^^° with is being the mean distance between particles. 
For screened by magnetic field spins C,o is to be assumed phenomenological and taken from 
experiment. 

We establish the explicit form of Eq. (ITOl) . The growth of boson spin fluctuations should 
cancel term, Eq. ([5]), in the Hamiltonian. Therefore, for dense {vo < C,o) Bose gas there 
should be s^+s = —hbk/M. 

The Hamiltonian of bosonized infinite anyon Coulomb gas with interaction of spins has 
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a form 
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For anyon Coulomb gas with density parameter > 2, where is To in Bohr radius as 
units, the approximate ground state energy [sl per particle 

Eo _ cl'cpl\v.rs) , lvi'l\v.rs) .... 



expressed in i??/, Rydberg energy units, and at u = 1 with 15% accuracy represents Quantum 
Monte Carlo result Q of Tanatar and Ceperley for 2D spin polarized electrons. Here 



f{u, r.) ^ u^^l - r.)e~- + ^ ^"""r^ • (13) 
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In Eqs. ([T2D and ([T3]) cvkc = 3.2903 and cbl = 1.2934. 

The expectation value calculation of energy with Hamiltonian, Eq. (fTTl) . gives analogous 
Eq. ( fT2i) expression for gas of bosonized anyons, however, one needs to replace in it u by 
j/g-'f-s/^o and now is in units. 

III. COMPARISON WITH EXPERIMENTAL PHASE DIAGRAMS 

Microscopical mechanism of High-Tc superconductivity (HTSC), as also the symmetry of 
its condensate wave function, are challenging problems of condensed matter physics. Some 
aspects of the theory, related to these topics, have been considered in papers [l^. We 
apply the model for clarification of experimental phase diagrams of hole and electron doped 
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cuprates proposed recently in Ref. [ll| and Refs. 

In our treatment for the description of HTCS we consider the bosonized fermions with 
u = 1. To become fermions bosons should overcome the energy difference 

,B 7(l-e--/^°)r/^(0,r.) 

o 2/3 8/3 ' y ) 

the gap of superconductivity. Our approach, as in jsl, corresponds to spinless or fully spin 
polarized fermions. One needs to have deal with normal, i.e., no spin polarized electron 



5 



liquid. However, the accuracy of our calculations is lower than the difference of Tanatar and 
Ceperley data for ground state energy for these both phases of electrons. 

We express as function of density of dopants rig = l/(7rr^) and connect with p, 
fractional part of doped hole or electron per atom Cu. At big values of or small Ug one can 
neglect the exponential factor in Eq. fll4l) and Hs ^ p. At small or big n^, without 

this factor, Af would have Af ~ n^J"^ ~ p^/^, but e"**^/^" suppresses this dependence to zero 
and the law of it depends from function S,o{p)- For this limit of we assume that A^(p) 
coincides with experimental dependence Eg{p). Extrapolating this asymptotic expression 
of Af (p) to small values of p and equating it to Eg{p) one finds the empirical dependence 

Up)- For it Up) -'^/Egip). 

To be sure in the correctness of our 2D density of holes, we ex pre ss it in cm~^ units and 
compare with experimental one. For elementary structural cell [ij] of almost all cuprate 
ab planes a = 3.8lA and b = 3.89A. Assuming that it has one atom of Cu, the density 
is Uab = Nab ■ P cm~'^, where Nab is number of elementary cells per Icm^ square. One 
obtains Uab = 6.7472 ■ p ■ 10^^ cm~^. We find experimental density n'^^ from paper [l^ 
for Y-123 (YBa2Cu30i^s) compound at optimal doping concentration of holes p ~ 0.16. 



It is ra^^^ = 0.9137 ■ 10^^ cm ^. Our optimal doping value for Uab is 1.0795 • 10^^ cm 
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From experiment [15l| also leads important information about approximate equality of total 
density of holes to density of HTCS carriers. A comparison of this value for Uab with 
experimental concentration 16| of carriers for Fractional Quantum Hall Effect (FQHE) 
npQHE ~ 10^^ -i-10^^ cm~^ shows importance for HTCS, as also for FQHE |6;], of long-range, 
not screened. Coulomb potential interaction between particles. The screened potential is 
supposed to be as justification [17] for the treatments based on the Hubbard model. 

The Fig. 1 displays pseudogap (PC) boundary energy Eg (Fig. 11 from paper 11|). 
HTCS gap energy Ao = 4/^^ T^, which was evaluated by empirical formula = Tc^maxi^ — 
82.6(p - 0.16)2] with Tc^rnax = 95 K for Bi -2212 {Bi2Sr2CaCu20s+5) compound, and 
HTCS gap energy calculated from Eq. (|T4l) as function of p. As we see our A^ has the 
same magnitude as experimental gap, but is qualitatively different from generally accepted 
"dome" like temperature-concentration phase diagram. However, it is in accordance with 
Fig. 10 of paper ll(| of Tallon and Loram and their conclusion that PC energy Eg up to 
Pc ~ 0.19 separates Bose-Einstein condensate into regions, where density of Cooper pairs is 
small and big (weak and strong superconductivity). 
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FIG. 1: The experimental PG Eg, HTCS gap Aq = AKb Tc (experiment for hole doped Bi - 2212 
compound), and calculated from formula Eq. (|14p HTCS gap for bosons energies in Kelvin 
temperature (K) units as function of concentration of holes p. 



For phase diagram data of electron doped cuprates we use Ref. 



12| for NCCO 



{Nd2-xCexCuO^ and Ref. [13| for PCCO {Pr2^xCexCu04). It was shown that 
Eg/{kBT*) ^ 10 for NCCO in 0.05 < j9 < 0.10 and Eg/{kBT*) ^ 11 for PCCO in 
0.11 < p < 0.17, therefore, we assume Eg/{kBT*) ^ 10 for entire interval of p and in- 
terpolate Eg{p) with dependence Eg = — 3916.67p + 675.83. For experimental HTCS gap 
we also assume Ao/{kBTc) ~ 10. Fig. 2 shows the p dependence of experimental Eg, 
Ao = 10k bTc and Af calculated from Eq. (IT^ by using the above spacing constants of a 
and b for elementary structural cell. Comparing with Fig. 1, we see the same qualitative 
and quantitative result. More obvious is extension of our Af to small values of p, while 
experimental A^ starts with p = 0.13. However, absolute values of both HTCS gaps are 
nearly equal. 

The spin correlation length sharply increases when p approaches pc, which might 
mean the vicinity of phase transition (PT). Fig. 6 of Tallon and Loram paper shows that 
experimental short-range Anti- Ferromagnetic (AF) correlations scale like Eg{p) dependence 
and vanish at p^. The presumption about correlations of spins of bosons would mean that 
there would exist the competing of these correlations with AF correlations when increasing 
of correlation length C,o leads to revealing the Fermi like spin correlations of bosons in the 
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FIG. 2: The experimental PG Eg, HTCS gap Aq = WKb Tc (experiments for electron doped 
NCCO and PCCO compomids), and calculated from formula Eq. (|14p HTCS gap for bosons 
energies in Kelvin temperature (K) units as function of concentration of electrons p. 

extended region of sample and thus, assuming that these spins now interact with ones of 
fermions, to suppressing of short-range AF correlations inside of this region. The definition 
of competing of HTCS and AF phases, which is widely used in the literature, in this case 
would find the natural interpretation. Only two phases, which have the relative nature, can 
compete with each other. AF magnetic phase and based on Cooper pairs HTCS phase can 
not coexist, at least, for high, as for HTCS, temperatures, because they are antagonistic. 
One would expect the Cooper pairing in right side from p^, where the magnetic correlations 
are not strong as AF correlations. 



We do the hypothesis that due to vicinity 14| of structural PT to HTCS state in the 
cuprates for concentration of holes below optimal the induced by it mechanical strain would 
strengthen a quadratic striction and therefore, the PT into Bose-Einstein condensate is not 
a second order, as it should be, but first order, close to second one The PG regime 

would correspond to meta stable phase of bosons. At boundary Eg bosons finally undergo 
PT into fermions. The critical concentration pc, at which Eg = , would be considered as 
critical point of first order PT. In general, the gas of doped holes and electrons in the PG 
region might be considered as mixture of Coulomb interacting single particle bosons and 
fermions. In the deep underdoped regime and low temperatures dominate bosons, close to 
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Eg dominate fermions. 



IV. SUMMARY 

We have shown that interaction of spins Sz — h/2 of anyons or fermions with statistical 
magnetic field leads to bosonization. For fermions, as the particular case of anyons, we have 
applied our model for the possible clarification of phase diagrams of HTCS of hole and elec- 
tron doped cuprates for region below preudogap Eg and found qualitative and quantitative 
agreement with experiment. 

B. A. gratefully acknowledges the Korean Research Foundation Grant (KRF-2004-005- 
C00044) for support of work. 
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